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The equations describing isothermal diffusion in the dilute range of the ternary system 
Ca(OH) 2 — H3PO4 — H 2 are derived in two ways, first, assuming that the components are electro- 
neutral species and, second, considering the actual ionic species present in solution. It is shown that 
the two models are thermodynamically equivalent. The theory permits the calculation of the four 
fundamental diffusion coefficients (phenomenological coefficients) in the concentration range where 
the Debye-Hiickel theory suffices for the calculation of ionic activity coefficients. The equations can 
be used to test the Onsager reciprocal relations for the diffusion process in the above system. The 
ionic model was used to calculate practical diffusion coefficients for the electroneutral components 
from the limiting equivalent conductances of the ions in solutions saturated with respect to hydroxy- 
apatite, Caio(OH) 2 (P04)<;. Large diffusion interferences, as revealed by relatively large values for the 
cross-terms in D;j, are predicted even for solutions with total molarity in the order of 10~ 5 . Therefore, 
diffusion models based on independent fluxes of the components appear to be invalid. 
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1. Introduction 

Diffusion processes in dilute solutions of the ternary 
system Ca(OH) 2 — H. { P04 — H 2 have biological and 
theoretical importance. Mechanisms for calcification 
and for formation of caries in teeth based upon 
diffusion of calcium and phosphate species have been 
advanced [1, 2, 3]. 1 Ternary systems are the simplest 
systems in which the validity of the Onsager's recip- 
rocal relations [4, 5, 6] can be tested for a diffusion 
process. The theoretical and experimental treatments 
of the present system have, besides the difficulties 
inherent to any multicomponent system, the complica- 
tions originating from the presence of a weak electro- 
lyte (phosphoric acid) and the impossibility of using 
optical methods to monitor the diffusion process such 
as has been done at higher concentrations [7, 8, 9]. In 
order to test Onsager's reciprocal relations, it is neces- 
sary to calculate the partial derivatives of activity 
coefficients with respect to the concentration of the 
components. Previous treatments of this problem in 
ternary systems have either assumed invariance of 
activity coefficients with concentration [10] or have 
used assumptions applicable to symmetric electro- 
lytes [11]. The present treatment applies to the 
ternary system Ca(OH) 2 — H 3 P0 4 — H l O in the 
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concentration range where the Debye-Hiickel theory 
has been shown [12, 13] to describe adequately ionic 
activity coefficients. Derivations will be given for two 
models in which electroneutral components and ionic 
constituents are considered, respectively. It will be 
shown that these two models are thermodynamically 
equivalent. 

The fnathematical treatment is applied here to the 
calculation of "fundamental" and "practical" diffusion 
coefficients for electroneutral components in solutions 
saturated with respect to hydroxyapatite; such sys- 
tems were selected because of their biological rele- 
vance. In a subsequent publication, the present 
derivations will be applied to results obtained in diffu- 
sion experiments using diaphragm cells. 

2. Procedure 

2.1. Selection of Model 

All the equilibrium properties of the ternary system 
referred to in this paper can be defined in terms of 
the components Ca(OH) 2 , H3PO4, and HjD. For 
diffusion purposes, the number of components, C", 
may be obtained from the expression C = n — r, in 
which n is the total number of chemical species, ionic 
and electroneutral, present in the system and r is the 
number of restrictions operating on those species. 
Usually, for the case of strong electrolytes the only 
restriction is electroneutrality. If Ca(OH) 2 , H,iP04, and 
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H 2 are selected as components, in general, n would 
comprise eight constituents: H 2 0, H3PO4, H 2 POj, 
HPOj, POf , H + , Ca+ + , and OH". The system would 
have five restrictions: the three expressions for the 
three ionization constants of H3PO4, the expression 
for the ionization constant of H 2 0, and the condition 
of electroneutrality. Therefore, the system is ternary. 2 
It is neither fruitful nor realistic to consider all the 
foregoing consitutents and restrictions. The derivations 
would become unnecessarily complicated; besides, 
because of the magnitudes of the ionization constants 
of H3PO4, in any given system only two phosphate 
species can be present in significant concentrations. 
Furthermore, as done by other investigators [10, 14], 
it will be assumed that the solvent, H 2 0, is undisso- 
ciated. If now, neutral components that do not contain 
the ion OH - are selected, n will be equal to five 
(H 2 0, two phosphate-containing species, Ca ++ , and H + ), 
and r will be two (equilibrium between the two 
phosphate-containing species and electroneutrality); 
hence the simplified model is ternary for diffusion 
purposes. The two phosphate species considered here 
are H 2 P04~ and HPOf which are the only species 
present in significant concentrations in dilute solutions 
and in the pH range 4.5 to 7.5. The other two com- 
ponents, besides H 2 0, correspond to electroneutral 
combinations of the ions considered and their selection 
is arbitrary. For convenience, CaHPCX and H3PO4 
are selected here, but the treatment would apply 
equally well to, for example, CaHP0 4 and Ca(H9P04) 2 . 

2.2 Derivations 



For the foregoing reasons, it is assumed that the 
derivations shown here, although strictly applicable 
only to a solvent-fixed frame of reference, are valid 
for the systems under consideration, since their high- 
est concentration is in the order of 10~ 2 M. 

a. Electroneutral Model 

According to the principles of nonequilibrium 
thermodynamics [17, 18], the diffusion flow for the 
component i, J/, in a ternary system may be described 
by 



Ji = -^L ij {dfi J ldx) 1 = 0,1,2 



(1) 



in which L\j are phenomenological coefficients and 
the quantity within parenthesis is the gradient of the 
chemical potential fXj for component j along x; the 
summation in (1) includes the solvent designated by 
the subscript zero. 

From the Gibbs-Duhem equation, at constant tem- 
perature and pressure, it is obtained 






(2) 



Solution of eq (2) for 5 [lo/dx and substitution into eqs (1) 
gives 

J i = -j^L ij (dtJL j /dx) i=l, 2 (3) 



For the range of dilute solutions considered here, 
total molarity in the order of 10 4 to 10 ~ 2 , it is reason- 
able to assume that the partial molar volumes of the 
components are independent of concentration. Under 
these conditions, the practical diffusion coefficients, 
(Dij) c , obtained in a cell-fixed frame of reference (a 
porous diaphragm cell), become identical with those 
referred to a constant- volume frame of reference, 
(Djj) v [10, 15]. Furthermore, the solvent-fixed diffusion 
coefficients, (Aj)o, are related to the volume-fixed 
frame of reference [10, 14] by 

i= 1 ,...., q k = 1 ,...., q 



in which L ij = Lj j — Cj/C ( )Ln ) . In eqs (3), the flow of the 
component i is given in terms of the chemical potential 
gradients of the solutes; therefore, the phenomeno- 
logical coefficients, usually called "fundamental 
diffusion coefficients" £#, display Onsager's reciprocal 
relations [14]. 

The flow Jt in eqs (3) may now be given in terms of 
practical diffusion coefficients, Djj. In a ternary system, 
only the concentrations, d and C 2 , of two of the 
components are independent variables. Therefore, 



dx k~ t dCk dx 



7=1,2. 



Substitution of eq (4) into eq (3) gives 



(4) 



in which d and C are the concentrations of the ith 
components and the solvent, respectively; Vq and Vk, 
the partial molar volumes of the solvent and the com- 
ponent K, respectively; the summation is taken over q 
components, excluding the solvent. It is apparent that 
at infinite dilution, (Dy)o= (Dij) v ; even for concen- 
trations about 1M, the two sets of coefficients do not 
differ by more than 1 to 5 percent [16J. 



2 In this paper, the presence of calcium-phosphate ion pairs is ignored; in the concentra- 
tion and pH ranges for which the present derivations are intended, the concentration of 
these ion pairs is small. Introduction of these species, however, does not alter the ternary 
nature of the system since for each one there is a corresponding dissociation expression 
(restriction). 



, __4> j 4* dnjdCk 
Jl ^ Lii AdC k dx 



J=l 



or 



''--?. D »f 



i=l,2 



(5) 



(6) 



in which the practical diffusion coefficients, D,a are 
given by 



D *-l L » dCk 



i = l,2 

£=1,2. 



(7) 



56 



The four practical diffusion coefficients D^ can be 
obtained through a minimum of two diffusion experi- 
ments by the use of eqs (6) with the appropriate 
boundary conditions. For the case of diaphragm cell 
experiments, however, the integrated form [19] 
of eq (6) does not permit an analytical solution for the 
Dji : ; the calculation of these coefficients in ternary 
systems has been effected by either imposing restric- 
tions in the experimental conditions [20, 21] or approxi- 
mations in the flux equations [22]. Recently, statistical 
procedures that obviate these limitations have been 
reported [23, 24]. 

From eqs (.7), the following relations may be ob- 
tained by solving for the Lif. 



L n = (DnlUL22 — Dv2lJL2i)ls 

L 12 = (D V1 fJLxx ~ D\ xfJL i2 ) /s 
L 2 \= (D 2 ifJL 22 -D 22 fl 2] )/s 

L 22 = (D 22 fjL U —D 21 /jl v >)Is 



(8) 



in which /Ay stand for the partial derivatives d/JLi/dCj 
and the denominator s is given by 

Equation (8), previously used by other investigators 
[14, 25] permit the calculation of the fundamental 
diffusion coefficients, Ly, from practical diffusion 
coefficients that are determined experimentally, 
provided that the values of the various partial deriva- 
tives are known. The method for evaluation of these 
derivatives is given next. 

The chemical potential of the ith component, /jlj, is 
defined in terms of its activity «,, by 

/jlj = /jL°+RT\naj 

in which /x° is a constant under fixed temperature and 
pressure, and R is the gas constant. The partial 
derivatives, /JLij, are given by 



m 



RT 



Oij 



i=l, 2 
7=1,2 



(9) 



in which a t j represents the partial derivatives daj/dCj. 
Thus, the evaluation of fly is accomplished if appro- 
priate expressions of ay are available. 

As shown in eqs (46) and (47) in appendix 1, the 
activities, a x and a 2 , for components 1 and 2 (CaHP04 
and H3PO4, respectively) are functions of their con- 
centrations C\ and C>. Other quantities in those 
expressions (hydrogen activity and ionic activity 
coefficients) are also functions of C\ and C>. Therefore 



dij 






Mr£l 



dCj/ ,- k , II. V, 

y ( da <\ ( d y\ 

f\dyi)c j ,c k , ll \dCj) c 



dai 
dH 



0Cih k 



i = l,2 

7=1,2 
£=1,2 
j*k 



(10) 



in which H represents the ionic activity of hydrogen, 
and y/ represents the activity coefficient of each of the 
/ ions assumed to be present in the system (i.e., 
Ca ++ , // + ,H 2 P0 4 -, and HPOf). 

The electroneutrality function E, eq (42) in appendix 
1, which will be used later, is defined on the basis that 
only these ions are present. In the subsequent equa- 
tions through eq (19), indexes i, j, and /r, assume the 
same values as in eq (10). 

In the dilute range of the ternary system under 
consideration, where the Debye-Hiickel theory pro- 
vides an adequate expression for the calculation of 
ionic activity coefficients [12, 13], ji are explicit 
functions of the ionic strength/. Therefore, the second 
factor in the summation term of eq (10) becomes 



dyi 
dCjJ Ck 



dyi / dl 
' dl \dCj/ Ck - 



(ID 



By substitution of eq (11) into eq (10) and making 

dyi 



Y \dyiJcj,c k ,n dl 



eq (10) becomes 
_/ den 

aii 'Uc-j 






+ a? 



y t \dCj Jc 



dl 

dCj 



(12) 



By differentiation of eq (43) in appendix 1, it follows 
that 

dCj) Ck \dCj} Ck , », y,\dH) Cj , Ck , V/ \dCj} Ck 



+2 



dyi/c 



dyi 



(13) 



By substitution of eq (11) into (13) and making 

^\dyi/cj,c k ,HdI 

eq (13) becomes 

dCj/c k \dCj/c k ,H,y l 

-) 



+ 



(K) (*n\ +I *(BL) 

\dH) tu e t .y.\BC J L UCjJc.. 



(14) 



Solving eq (14) for (dJ/dCj) Cfc 

'*) 

'j/c k ,H,V l 



dl_ 



dl_ 



+ 



SI 



»L 



5/7 
dCj 



(1-1*)- 



(15) 
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Substitution of eq (15) into eq (12) gives an expression 
for aij in terms of partial derivatives that can be readily 
calculated from experimental measurements, except 
for the derivative (dH/dCj) Cf .. The derivation for 
(dHldCj)c k is given next. 

Differentiation of the electroneutrality function E, 
eq (42) in appendix 1, gives 

\dCje k \dC J )e k .».r l \dH)r J .c k .y,\dCj)c k . 

r\dWc J ,e k .n\iCjc k . (16) 
By substitution of eq (11) into eq (16) and making 

S3E\ dg =E * 

Y Kdyi/n.u.H dl 



eq (16) becomes 
/ BE \ „ (BE 



.dCj 



= = 



dCj 



\dH/c J ,c^y l \dCj) c 



+ E* 



dl 

dCj 



(17) 



Substitution of eq (15) into eq (17) gives the solution 
for (dHldCj) Ck . 



dH_ \ 

dCj )c 



(1-/*) 



dE 



+ E* 



dl 



dCj/ Ck ,ii.y l V dCj Jc k , ii, y, 



(18) 



b. Ionic Model 

In the foregoing treatment, the components have 
been considered as neutral species, although in defin- 
ing their activities, their ionic nature was taken into 
account. Consideration of an ionic model for the 
present system is advantageous in that practical diffu- 
sion coefficients can be calculated for the electro- 
neutral components on the basis of ionic conductances 
[10]. It is then possible to compare the practical diffu- 
sion coefficients derived here with those obtained 
experimentally. 

The ionic model selected here is given in eqs (20) 
to (22) 



iHP0 4 ^Ca ++ + HPOT 


(20) 


H: S P04^2H + +HP04= 


(21) 


H2PO4- ** H + + HPOf 


(22) 



This ionic model implies that there is no undissociated 
H3PO4, an assumption which is justifiable for the pH 
range considered here. It will be shown now that this 
system is thermodynamically equivalent to the one 
in which CaHP04 and H3PO4 are considered as the 
neutral components. In the following derivations, the 
numerals 1 and 2 are used for the constituents CaHPQi 
and H3PO4; the numerals 3, 4, 5, and 6 represent the 
ions H2PO4-, Ca + + , HP0 4 = and H + , respectively. 

The stoichiometric coefficients in the right-hand 
terms of eqs (20) to (22) are represented by Pij. For 
example, V\4 is the stoichiometric coefficient for the 
Ca ++ that comes from CaHP0 4 and v^ is for the 
HPOf that comes from the CaHP0 4 . The numerical 
values for the Vifs are: Vu = ^15 = ^25 = ^5 = ^36 — 1; 

*>16== ^24 = I>34 =0; V 2 6 = 2. 

On the basis of the ionic model, the entropy produc- 
tion, crT, in the diffusion system is given by [18, 26] 



Substitution of eq (18) into eq (15) gives 



dl 

dCj 



dl 



m 



dCj J c k , h, Yf V dH 1 cj, c k , y l 



) 



di 



dH J c c y 



ML 



*.H,y, 



(19) 



{i - n ^K^ 



\8H / Ci ,c L .,y, \BH )c jt c ki y l 



The right-hand terms in eqs (18) and (19) involve 
quantities that can be readily calculated from experi- 
mental measurements. The actual expressions used 
for such calculations are given in appendix 2. Substitu- 
tion of eqs (11), (18), and (19) into eq (10) gives the 
needed expression for the evaluation of a,,; then, by 
the use of eq (9), evaluation of /x,j is accomplished, and 
substitution of the fly into eqs (8) permits the cal- 
culation of the fundamental diffusion coefficients, 
Ltj, and, therefore, the testing of the Onsager's 
reciprocal relations. 



<tf=2JA = 2 /&+**» 



(23) 



j=3 



j = 4 



in which Jj represents the diffusion flow of theyth ion 
and Xj its generalized driving force (in this case, its 
electrochemical potential gradient); cr is the rate of 
entropy production at temperature T. 

We consider now as new generalized forces the 
chemical potential gradients, Xu of the neutral com- 
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ponents 1 and 2 and the electrochemical potential X :i . 
The transformation equations for the new forces are 



i=4 



(£ = 1,2) 



(24) 



which are linear combinations of the old forces Xj. 
The transformation equations for the new flows, derived 
from the mass balance condition in eqs (20) to (22), are 



2 

I 



Jj = 2 v vJ* ~ v vJz 0' = 4 to 6) . 



Substitution of eq (25) into (23) gives 



(25) 



1=1 j = 4 ^j=4 / 

But, according to eq^ (24) the summation term within 
parenthesis equals A3, therefore, the quantity within 
parenthesis vanishes and 



crT=2jiXi 



(27) 



hence, the ionic model is thermodynamically equivalent 
to the diffusion system in which the neutral compounds 
CaHP04 and H3PO4 are considered as components. 

The procedure to calculate fundamental diffusion 
coefficients for neutral components from ionic con- 
ductances has been given in detail by Wendt [10] and 
it was tested by the same author with diffusion data [27] 
for the system H2SO4 — Na 2 S04 — H 2 0. For this reason, 
only a summary of our derivation is given next in which 
we apply Wendt's treatment to the present system. 
We use the same terminology except that his Ly 
(for ions) corresponds to our Ty and his a X) correspond 
to our (Xij. 

The diffusion flows Jj of the j ions are related to the 
chemical potential gradients, X H , of the two neutral 
components by 



In eq (29), 8/* is the Kronecker delta (it has the value 
1 for 7 = /c, and zero for 7 4= A), and the Ps are the main 
terms of the fundamental diffusion coefficients for 
the subscripted ions, 3 to 6. Equation (29) is consistent 
with Wendt's assumption [10] that the coefficients Tjk 
are zero for j =1= k. These ionic coefficients can be 
calculated from 



i> 



(kjCj)l I Zj I F 2 X 10 7 y = 3 to 6 (30) 



in which Xj is the limiting equivalent conductance of 
the ion 7, Cj its molarity, | Zj | the absolute value of its 
valence and Fis the Faraday constant (96,493 coulombs 
per equivalent of charge). Values used for Xj in this 
investigation are given in appendix 1. 

Expanding eq (25) and solving for J\ and J> (using the 
numerical values for the corresponding i/y) we get 



Ji=A 

/t=i(/s+/e). 



(31) 



The expressions for the ionic flows J :i , J a, and J% 
given by eq (28) can now be substituted into eqs (31) 
to obtain 



ii — 2 2 OLkrXr&tk 



(32) 



j * = \ i i "*x&k+l i 2 otoX&k (33) 

ti = 1 A- = .1 ^ » = 1 ie = a 



n=\ A = 3 " n = \ A=3 

But from eq (3) we can also write 



J\ / J L\ n X„ 

n= 1 

2 

72 — — 2j L>2nX n > 



(34) 
(35) 



Through a comparison of the coefficients of X n in 
eqs (34) and (35) with those in eqs (32) and (33), respec- 
tively, it is obtained 



Lln = ~ 2 0LknQ>4k 



A-r'2 2«taWV 0=3 to 6) (28) ^ 1 A ^j^ + ft*) 

n=\ A=3 /. *-** 



n=\,2 (36) 



in which v is the value of the determinant of the v\ 3 
matrix (equal to — 1 in the present case); «a-w, the 
elements of the matrix [a* n L have the numerical 
values a :u = a>,y> = ol-,\ = a^i = 0, CH33 = a 4 i — 0L42 — ol%2 = ~ 1 , 
a->2 — ot%z— 1, «43 = 2, tt53=— 2; the coefficients Cljk 
are given by 



tyk — 5; a- IV — 



6 
g*3 



7 = 3 to 6 
k = 3 to 6. 



(29) 



Equations (36) together with (29) and (30), permit the 
calculation of theoretical fundamental diffusion co- 
efficients, Lij 9 from a knowledge of the limiting con- 
ductance of the ions involved. 

3. Calculation of Diffusion Coefficients 

The values for Ly and Aj, calculated from eqs (36) 
and (7), respectively, are shown in table 1, together 
with the composition of the 13 solutions used for 
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this purpose. The concentrations C\ and C> were 
taken from the solubility isotherm for hydroxyapatite 
[13]. Since, in the present derivations, ion pairs are 
ignored, the values for pH in table 1 were calculated 
solving eq (42) for (H + ). The calculated pH values 
do not depart appreciably from the corresponding 
values for the isotherm [13]. 

The data in table 1 suggest that, while the inter- 
ference represented by D vl is important in all the sys- 
tems considered, the effect of Ai becomes important 
only in the systems having pH's lower than about 6. 
Study of system number 4 showed this to be the case. 
Calculations of the/)/, were made keeping the concen- 
tration of one of the components constant and decreas- 
ing the concentration of the other component. The 
results of these calculations are shown in figures la 
and lb. It is apparent that, as the concentration of 
CaHP0 4 , Ci, approaches zero, figure la,Z) 2 2 approaches 
a value very close to 1.575 X 10 -5 ; this is the Nernst 
limiting diffusion coefficient calculated with the values 
for the limiting equivalent conductances given in 
appendix 1 for the ions H + and H 2 POf; in other 
words, the diffusion of 4.888 X 10" 5 M H 3 P0 4 in the 
binary system H3PO4 — H 2 essentially fulfills the 
expectations for infinite dilution. 

The calculated Nernst limiting diffusion coefficient 
for CaHP0 4 is 1.342 X 10 5 . Examination of figure lb 
indicates that this limiting value is not approached 
even at the lowest concentration of H3PO4. It must 
be concluded, therefore, that, according to the model 
adopted here, the interference of H3PO4 on the dif- 
fusion of CaHPCX persists even at very low concen- 
trations of both components. This conclusion was also 
verified by similar calculations on systems more dilute 
than system 4. 

The values in table 1 and the curve in figures la and 
lb clearly show that in diffusion processes occurring in 
dilute solutions of the system CaHP0 4 — H3PO4— H 2 0, 
the magnitude of the interferences, as given by the 
cross terms in Dfy, will affect the fluxes appreciably and 
therefore cannot be ignored. Also, for a test of the 
Onsager reciprocal relations, experimental conditions 
should be selected such that the interferences are 




lb 

C 1 -9092x10" 5 «i/i 



C.xlO niol/l 



FIGURE la. Plot of the main diffusion coefficient 0/H3PO4 against 
concentration o/CaHPO.*. 

The concentration of H 3 P0 4 is fixed at 4.888 X 10 5 mol/1. 

FIGURE lb. Plot of the main diffusion coefficient of CaHPO 4 against 
concentration of H3PO4. 

The concentration of CaHP0 4 is fixed at 9.092 X 10 ~ 5 mol/1. 



Table 1. Calculated Ui and D u /or the Ternary System CaHP0 4 (l)—H :i P0 4 (2)—H 2 at 25° 



Systems 


Ci, mol/1 X10 4 


C 2 , mol/1 X10 4 


pH 


L n xi0 i:i 


L 12 xiO ,:5 


L n x 10 13 


L, 2 X 10 13 


D n x 10« 


D 12 xi0 fi 


D, x X 10« 


D, 2 x 10" 


1 


0.319 


0.0596 


7.521 


0.4201 


0.05890 


0.05890 


0.1002 


6.552 


1.509 


0.05540 


8.724 


2 


0.420 


0.114 


7.311 


0.5468 


0.1139 


0.1139 


0.1851 


6.533 


1.536 


0.08048 


8.644 


3 


0.658 


0.280 


7.008 


0.8343 


0.2863 


0.2863 


0.4288 


6.506 


1.564 


0.1085 


8.533 


4 


0.909 


0.489 


6.808 


1.131 


0.5079 


0.5079 


0.7160 


6.493 


1.568 


0.09988 


8.483 


5 


1.54 


1.07 


6.512 


1.857 


1.129 


1.129 


1.458 


6.500 


1.536 


-0.008148 


8.497 


6 


2.71 


2.20 


6.220 


3.198 


2.370 


2.370 


2.846 


6.568 


1.423 


-0.3150 


8.704 


7 


4.03 


3.51 


6.027 


4.703 


3.796 


3.796 


4.401 


6.674 


1.275 


-0.6894 


9.008 


8 


6.74 


6.20 


5.785 


7.790 


6.752 


6.753 


7.585 


6.920 


0.959 


-1.448 


9.670 


9 


10.3 


9.75 


5.598 


11.85 


10.65 


10.65 


11.76 


7.217 


0.5940 


-2.303 


10.44 


10 


19.9 


19.3 


5.312 


22.79 


21.19 


21.19 


23.01 


7.982 


-0.2990 


-4.300 


12.28 


11 


31.4 


30.8 


5.113 


35.86 


33.81 


33.81 


36.46 


8.851 


-1.270 


-6.385 


14.26 


12 


64.0 


63.5 


4.787 


72.96 


69.62 


69.62 


74.83 


11.23 


-3.822 


-11.64 


19.32 


13 


105. 


105. 


4.485 


120.2 


114.6 


114.6 


124.6 


15.78 


-8.510 


-20.87 


28.42 



Units of L u : Mol 2 /(J • cm • 
Units of Djj: cm 2 /s 



s) 
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detectable with the measurements used within 
reasonable errors. For example, with system 1 of 
table 1 it would be practically impossible to arrive at 
a reliable value for L 2i . On the other hand, it is pre- 
dicted that diffusion experiments with systems 12 or 13 
will yield the required information for calculation of Ly. 

4. Appendix 1 

C\ = molar concentration of CaHP0 4 

Ci = molar concentration of H 3 P0 4 used in prepara- 
tion of solution 

P\ = molar concentration of H 2 PO,r 

P> = molar concentration of HPOf 

yi = molar activity coefficient for H2PO4" 

y 2 — molar activity coefficient for HPOr 

y-.] = molar activity coefficient for Ca ++ 

y 4 = molar activity coefficient for H + 
[H] — molar concentration of H + 

H = activity of H + 

K] = first ionization constant of H t} P0 4 , 7.10 S X 10~ 3 
[28] 

K> = second ionization constant of H, } P0 4 , 
6.33 H xl0- 8 [29] 
(P 2 ) = activity of HPOf 

A.3= limiting equivalent conductance of H 2 P0 4 ", 32.3 

k°4— limiting equivalent conductance of Ca ++ , 59.50 
X.5 = limiting equivalent conductance of HPOr, 43.7 

1 32 ! . 

k°6= limiting equivalent conductance of H + , 349. 81 
[31]. 

The total molar concentration of phosphorus in the 
solutions of CaHP0 4 in dilute H ;{ P0 4 is given by 
(C1+C2). The balance equation for phosphorus, for 
solutions with pH values below 7.5, may be written 



(C,+C«)=/>i+IY 



(37) 



The quantity P\ may be written in terms of P 2 by the 
use of the expression for the second ionization of 
H3PO4 ; then, solving eq (37) for P 2 gives 



Pi=(Ci+C,)/t* (-+#-). 

Vy2 K 2 yi) 



(38) 



Calling / 2 the quantity within parenthesis in the 
denominator of eq (38) 



f»,= (Ci+C,)/7*/i 



(39) 



and by substitution of eq (39) into the expression for 
the second ionization of H.-iPO-t, 

P 1 ^(C 1 + C i )Hlf 2 y 1 K 2 . (40) 

The function E, that defines electroneutrality, is 

H 



£ = = 2C,+ 2P,-Pi. 

7-i 



(41) 



Substitution of eqs (39) and (40) into (41) gives 



The ionic strength, /, is given by 

T-9C 1 H 1 C i+ C * 
7 " 2Cl + 2y 4 + 2/ 2 






H 4^ 
yiK 2 72. 



(42) 



(43) 



For brevity in the expressions given in appendix 2, 
we define the quantities within square brackets in 
eqs (42) and (43) as 



K 2 y\ y 2 



c=^- + - 4 



K>y, y,' 



(44) 



(45) 



The activities of CaHPOj and H :i PO t in solution 
are given by eqs (46) and (47), respectively. 



a x 



_ (a + C>C l )y, 

h 

_ (C t +Cj)H 



The function f t in eq (47) is defined by 

Hyt 7, 



(46) 

(47) 

(48) 



The ionic activity coefficients, y/, are calculated 
from the Debye-Hiickel theory 



In yr- 



AZfVI 
1+BaiVI 



(49) 



in which the constants^ and#, dependent on tempera- 
ture and the dielectric constant of water, have nu- 
merical values of 0.5092 and 0.3286 X 10 H , respectively 
[33] at 25 °C; Z\ is the valence of the /th ion, and a.\ its 
distance of closest approach. Numerical values for 
a x [34] are: Ca + + , 6X10~ 8 cm; H + , 9 X 10~ 8 cm; 
HPOf and H 2 P0 4 , 4 X 10« cm. 



5. Appendix 2 

In this appendix are given the expressions for the 
calculation of the partial derivatives used in the text; 
these expressions are based on the equations in 
appendix 1. 



dyi 
dl~ 



'-yi 



<*>(/) 



2I(l + BaiVT) 
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yhich 



$(/)=- 



di 



AZfyJl 

l+.Ba,VT 



= 2 + -^ 

dC\/ C 2,H,V, 2/~2 



dl\ = G_ 



ill 



l , (Ci+O) 



BH)c 1 ,c 2 ,y l 2y 4 2/2K2 



1 G_ 



a/ 



c, + c 2 



dy-Jc^H.v,^., 2f 2 y% 1/2 



-4 



a/_\ _ (Ci+c 2 ) 

Byi)c it c t ,H,y lfki 2/2^7? 



&-> 



dyJcuct.H.y^ 2y| 
daA _ (2Ci + C 2 )y 3 



3 C> /<•,,//. 7, /2 



3ai\ (C, 2 + CiC 2 )y3 /2ff J. 

a//J„, C2; y, _ /#tf 2 U. + y. 



/dai 



(C?+GC,)7iff 



Vayi/ccH.r^, /f/ky 



/da. 



(C?+C,C 2 )y 3 



\dy-2/r t .c 2 ,H,y l + . 2 Ziyi; 



daA = (C, 2 + C,C 2 ) 

Byz)c v C 2 ,H,Yj +3 /2 

dCo/c^ii^i \dCt/c 2 , H,y, fiKi 



daA _ (C, + C 2 ) / K 2 // 

dH) tvH . yi f?K> V' Hy 2 K t 



daA = (Ci+C 2 )H 

Byi)e v e v y l ^. 1 ftKiy? 



da*\ _ (C, + C i )K-> 

3yJo 1>% ,y J+2 y?K,y| 



dC\ Jc 2 ,H,yi ft 

dC2/c v H,y l h 



(—)■ 

\dH/e v c 2 ,y l T4 / 2 K 2 yi V / 2 



( ft + C,) ^_£ 



d£\ 



(C, + C2)// / 1 <? 



dyi)e v c i ,H,y l ^ fiKoy? \ f 2 

dE\ = (C + C 2 ) / 2 g 

372/0,, c 2 ,tf,y (+2 / 2 y| V /a 






c,,c 2 ,H,7 / + 4 



yf 
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